Observed Properties of Dark Matter: dynamical studies of dSph galaxies 
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The Milky Way satellite dwarf spheroidal (dSph) galaxies are the smallest dark matter dominated 
systems in the universe. We have underway dynamical studies of the dSph to quantify the shortest 
scale lengths on which Dark Matter is distributed, the range of Dark Matter central densities, and 
the density profile(s) of DM on small scales. Current results suggest some surprises: the central 
DM density profile is typically cored, not cusped, with scale sizes never less than a few hundred pc; 
the central densities are typically 10 — 20GeV/cc; no galaxy is found with a dark mass halo less 
massive than ~ 5.1O 7 Af . We are discovering many more dSphs, which we are analysing to test the 
generality of these results. 



I. INTRODUCTION AND METHODOLOGY 

Determining the nature of the dark matter is one of 
the key goals of contemporary astronomy. The extremely 
large mass to light ratios for certain Local Group dwarf 
spheroidal (dSph) galaxies suggest that these are the 
most dark matter dominated stellar systems known in 
the Universe. Given the apparent absence of dark mat- 
ter in globular clusters (length scale ~ 10pc), and the 
direct evidence that there is no dark matter associated 
with the galactic thin disk (length scale ~ lOOpc; [3,0]) 
the dSphs (characteristic radii ~ 300+ pc) are both the 
smallest systems to contain dynamically significant quan- 
tities of dark matter, and the most nearby systems where 
we may look for characteristic properties of the distribu- 
tion of dark matter (maximum density, etc ...), which 
may provide some knowledge of its nature. Addition- 
ally, simulations of dark halo formation in a cosmologi- 
cal context predict that dSph galaxies are (unmodified?) 
survivors from the earliest structures formed, so that un- 
derstanding their nature and structure also has general 
implications for galaxy formation and evolution. 

Is it possible reliably to determine the dark matter dis- 
tribution in a dSph galaxy? Fortunately, some dSphs are 
astrophysically simple, are old, are apparently in equilib- 
rium, have no dynamically-significant gas, but contain a 
(small) number of stars, which provide ideal kinematic 
collisionless tracer particles. With the new generation of 
wide-field multi-object spectrographs on 4-8m telescopes, 
it is now viable to obtain sufficient high-quality kinematic 
data to determine the gravitational potentials in all the 
Galactic satellites. There are two types of analysis fea- 
sible with current datasets: a straighforward pressure- 
gravity balance analysis, based on the velocity dispersion 
moment, using the Jeans' equations, and an analysis of 
the full projected kinematic distribution function, using 
the full information set but requiring more sophisticated 
analyses. These analyses are made more robust, and sim- 
plified, by stellar population studies of the dSph galax- 
ies. Ultra-deep HST imaging for Draco and UMi |l7j 
show their stellar populations are indeed indistinguish- 
able from those of old metal-poor globular star clusters, 



and hence have M/Ly ~ 2. 

The Jeans' equations for a spherical stellar system lead 
directly to a model-independent mass estimator (see pj 
Eq. 4-55 & 4-56) requiring only knowledge of the velocity 
anisotropy and the true radial velocity dispersion (v 2 ). 
Assuming spherical symmetry, it is straightforward to 
obtain (v 2 ) from the line of sight velocity dispersion (v 2 ) 
using Abel integrals. In applying the Jeans' equation 
analyses, the spatially-binned dispersion profile is fit by 
an appropriate smooth function (figure 1), while the light 
is fit by a Plummer law, as this is an excellent fit to 
available star count data in all available cases (cf Q). 
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FIG. 1: Draco surface brightness profile, with a fitted Plum- 
mer model, and the observed line of sight velocity dispersion 
profile, from Q. 
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The Jeans' equation moment analysis involves two ra- 
dial functions, the mass density and a possibly vari- 
able anisotropy in the stellar orbits, this second gener- 
ating an anisotropic stress tensor. It is this anisotropic 
stress which gives the dSph galaxies their shapes, as di- 
rect observational limits exclude any angular momentum 
support against gravity. Equivalently, radially variable 
anisotropy is degenerate with mass (figure 2) making 
any deductions on the inner mass profile being cored or 
cuspcd in general model-dependent. Further information 
is needed to break this degeneracy, and fortunately is 
sometimes available, as we discuss below. In general how- 
ever, full multi-component distribution-function models 
are required to use the information in the data to break 
this degeneracy. 



20 



Radial 

Step function 
Isotropic 
Tangential 



a 

b 




FIG. 2: An illustration of the possible effects of stellar orbital 
anistropy variations on kinematic modelling, from fl^l . 



II. RESULTS OF DYNAMICAL ANALYSES 

Are dSph haloes cusped or cored? Debate continues 
to rage about whether the cusped haloes always created 
in CDM simulations arc in conflict with observations of 
rotation curves of Low Suface Brightness galaxies. The 
gas-free nature of dSphs makes them kinematically clean 
systems in which to test theoretical predictions. Stellar 
velocities may also be used to place constraints on the 
steepness of any possible central cusp, whether due to a 
black hole, the intrinsic physical properties of the CDM 
(eg 01 )> or possibly even CDM as modified by a central 
black hole (eg [13). 

As a general result, in all cases with sufficient data we 
rule out (King model) mass-follows-light models. King 
models are not an adequate description of these galax- 
ies, all of which have high mass-to-light ratios (quoted in 



solar visual band units), and extended dark matter ha- 
los. Even in the inner regions mass does not follow light, 
while including outer data commonly we find a most 
likely global mass to light ratio which is very high, being 
for Draco ~ 440, 200 times greater than that for stars 
with a normal mass function (Fig. 1). The Draco halo 
models favoured by the data contain significant amounts 
of mass at large radii, leading to the observed flat to ris- 
ing velocity dispersion profiles at intermediate to large 
radii. 
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FIG. 3: Derived inner mass distributions from Jeans' eqn 
analyses for four dSph galaxies. Also shown is a predicted 
r _1 density profile. The modelling is reliable in each case 
out to radii of log (r)kpc~ 0.5. The unphysical behaviour at 
larger radii is explained in the text. The general similarity of 
the four inner mass profiles is striking. 



In some cases (eg Draco, cf Figure 1) a decrease in 
the disersion is apparent in the outermost data. This re- 
mains to be fully understood. The decrease is apparently 
steeper than Keplerian, so does not indicate the (first) 
detection of the outer limits of the galaxy. It cannot in- 
dicate tidal perturbations by the Milky Way galaxy, and 
in fact robustly excludes such complications, as tides heat 
and do not cool (eg [Tfljh It may indicate a combination 
of a radially-changing orbital anisotropy and a complex 
stellar kinematic distribution function (Figure 2). In any 
case, this illustrates that an adequate understanding of 
the very outer parts of the dSph is still lacking, so that we 
are unable as yet to determine outer mass distributions, 
or 'total' masses. The situation in the inner regions is 
however better defined and understood. 

Figure 3 summarises Jeans equation models for sev- 
eral of the dSph, with in each case the simplest possible 
assumptions (isotropic radially-constant velocity distri- 
bution). It is apparent that the models are invalid at 
large radii, where an unphysical oscillation in the mass 
profiles is evident. In the inner regions however the fit to 
the data is good. In each case, a core-like mass distribu- 
tion is preferred. As noted above however, it is possible, 
by adding a radially-variable stellar anisotropy (essen- 
tially an extra function) to the fit, to fit steeper cusp-like 
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FIG. 4: TOP: Result of search for kinematic sub-populations 
in UMi. Coutours are linearly spaced stellar isopleths; the 
second peak of UMi's stellar population is visible at RA = 
12', Dec = 8' relative to UMi's center. Gray stars are UMi 
RGB member stars with measured velocities. The dots rep- 
resent points where a model with a kinematically cold sub- 
population is at least 1000 times more likely than a model 
composed of a single 8.8 km s - Gaussian. BOTTOM: Simu- 
lation of an unbound clump in a dark matter halo. The halo 
has a density law p(r) oc (a 2 + r 2 )~ 1 ^ 2 ; the time t is in units 
of Gyr. Top panel: a clump in a cored halo with a — 0.85 
persists for a Hubble time because the potential is nearly har- 
monic. Bottom panel: a clump in a cusped potential (a — 0) 
disrupts in less than 1 Gyr. The histograms at the upper 
right corner of each snapshot show the distribution of total 
velocity v = (v 2 + v 2 + v 2 ) 1 ^ 2 ; tick marks are spaced 1 kms -2 
apart. The stars in a cored halo remain coherent in velocity 
as well as in position. From @. 



central mass distributions. 

Can one distinquish between shallow and steep den- 
sity profiles using other information? Fortunately, in one 
special case, that of UMi, one can. In UMi, an other- 
wise very simple system from an astrophysical perspec- 
tive, an extremely low velocity dispersion sub-structure 
exists (figure 4a). We explain this |5| as a star cluster, 
which has become gravitationally unbound (the normal 
eventual fate of every star cluster, as mass is lost through 
normal stellar evolution), and which now survives as a 
memory in phase space. Why does it survive in config- 
uration space? The group of stars have the same mean 



velocity as the systemic velocity of UMi, so they must 
orbit close to the plane of the sky, and hence through the 
central regions of UMi. As figure (4b) illustrates, this 
is possible only if the tidal forces from the UMi central 
mass gradient are weak. In fact, survival of this phase- 
space structure in configuration space requires that UMi 
has a cored mass profile. Similar results are provided by 
the survival of the globular cluster system in the Fornax 
dSph 0. 



A. Distribution Function models 

Standard cosmological simulations tend to overpredict 
the numbers of low-mass haloes surrounding Milky Way- 
type galaxies compared to observations. One possible 
solution (eg ^J) is that all dwarf galaxies have very ex- 
tended haloes, and are therefore very much more mas- 
sive than previously thought. According to the Stoehr 
etal analysis, the observed dispersion profile of Draco and 
Fornax can be obtained by embedding the stellar distri- 
butions in ACDM haloes of masses > 10 9 M Q , extending 
to several kpc. While our existing data on Draco are com- 
patible with this hypothesis, it is not known whether it 
holds in general. More reliable masses and mass profiles 
are therefore desirable. 

The principal objective of our dynamical modelling 
work to date has been to move beyond the traditional 
King models for dSphs, in order to break the principal 
degeneracy of the problem, namely that a large veloc- 
ity dispersion at large radii in a dSph can be explained 
either by the presence of large amounts of unseen mass 
at large radii, or by tangential anisotropy of the veloc- 
ity distribution. We have therefore developed a family 
of dynamical models (described in detail in [l^) which 
span the full range of possible halo models from mass- 
follows-light to extended haloes. The models also in- 
corporate anisotropy in the velocity distribution. Our 
current modelling has focussed on spherically symmet- 
ric dark haloes and stellar distributions. Our analysis 
of our Draco data, together with extensive Monte Carlo 
simulations, has demonstrated that data sets of several 
hundred radial velocities are sufficient to discriminate be- 
tween halo models and to break the degeneracy between 
mass and velocity anisotropy provided these data are at 
sufficiently large radii The outermost data also 

place constraints on any systemic rotation or tumbling 
which might contribute, along with anisotropic kinemat- 
ics, to the flattening seen in some dsph galaxies, and 
provide further direct tests that a system is unaffected 
by tides, so that an equilibrium dynamical analysis is 
appropriate. 

The results of spherical modelling supports the basic 
results of the Jeans' equation analysis, though of course 
use more information, so are more robust. We are cur- 
rently implementing new families of models to allow more 
robust distinction between cored and cusped mass distri- 
butions in the dSph. One new option allowed by extant 



4 



DF models is to test the viability of the alternative grav- 
ity theory, MOND. Wilkinson et al. have carried out 
such an analysis, concluding that even with MOND dark 
matter is required. Available data for the Fornax dSph 
can in fact be modelled quite well by MOND, assuming 
the MOND scale parameter to be ao — 2.1 x 10~ 8 cms~ 2 . 
For UMi and Draco, however, even MOND requires ei- 
ther a high M/L value (~ 20), or a very different length 
scale parameter (ao ~ 50 x 10 _8 cms -2 ). 

III. GENERAL RESULTS: SYSTEMATIC 
PROPERTIES OF DARK MATTER ON SMALL 
SCALES 
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FIG. 5: Mass to light ratios vs galaxy absolute V magnitude 
for some Local Group dSph galaxies. The solid curve shows 
the relation expected if all the dSph galaxies contain about 
4 x 10 7 solar masses of dark matter interior to their stellar 
distributions. 



The profiles in figure 3, derived by Jeans' equation 
analyses, illustrate two of our basic results. In every case, 
the simplest analysis favours cored mass distributions. 
While cusped mass distribution can usually be fit to the 
data, in at least one case, UMi, there is very strong direct 
evidence that a cusp model is inadequate to explain all 
the available information. The conservative assumption 
is therefore that all the mass profiles are indeed cored, 
and are significantly shallower than r _1 . 



Secondly, all the dSph we have analysed to date 
show very similar, and perhaps surprisingly low, central 
dark matter mass densities, with a maximum value of 
~ 5 x 10 8 M Q kpc~ 3 , equivalent to ~ 20GeV/cc. Inter- 
estingly, the rank ordering of the central densities is in 
inverse proportion to system total luminosity, with the 
least luminous galaxies being the most dense. This is of 
the opposite sign to some CDM predictions, though we 
note a yet further complication involving system mass 
below, since most CDM simulations predict mass rather 
than luminosity. 

We have not yet detected a reliable Kepler ian decline, 
or any evidence for tidal truncation, in any dSph galaxy 
dark halo studied so far. We therefore (still) have no to- 
tal mass determinations. Nonetheless, the total masses 
determined out to the limits of the stellar distribution 
show an interesting systematic effect. This is shown in 
Figure 5, a plot whose style is adopted from [ljj. Fig- 
ure 5 displays the relationship between the luminosity of 
the lowest-luminosity Local Group galaxies against the 
derived (logarithmic) system mass-to-light ratio. In this 
parameter space, the solid curve shows the relationship 
anticipated if every dSph galaxy contains the same mass 
of dark matter - in this case some ~ 4 — 5.10 7 Mq within 
the volume which contains its stellar population. The 
system total luminosity is measured directly. 

It is apparent from Figure 5 that there is remarkably 
little spread in mass apparent among the galaxies with 
absolute magnitute fainter than ~ — 11. This relation was 
considered until recently to be a minor curiosity, since it 
covered the dynamic range only from My ~ —13 to —9, 
a mere factor of forty or so in luminosity, and included 
only 8 galaxies. However, the recent analysis [y| of the 
newly-discovered extremely low luminosity dSph galaxy 
UMa has extended the validity of the relation by another 
two magnitudes, now a factor of ~ 200 in luminosity, 
and to total mass-to-light ratios in excess of 1000. Sev- 
eral other new very low luminosity dSph satellite galaxy 
candidates have been discovered in the last few months, 
while new studies of several known galaxies have recently 
been completed 7J. It will be very interesting to see if 
these dynamical studies strengthen or disprove this ap- 
parent trend. 
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